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It is currently believed that the Standard Model is an effective low energy theory which in principle 
may contain higher dimensional non-renormalizable operators. These operators may modify the 
standard model Higgs potential in many ways, one of which being the appearance of a second 
vacuum. For a wide range of parameters, this new vacuum becomes the true vacuum. It is then 
assumed that our universe is currently sitting in the false vacuum. Thus the usual second-order 
electroweak phase transition at early times will be followed by a second, first-order phase transition. 
In cosmology, a first-order phase transition is associated with the production of gravity waves. In 
this paper we present an analysis of the production of gravitational waves during such a second 
electroweak phase transition. We find that, for one certain range of parameters, the stochastic 
background of gravitational waves generated by bubble nucleation and collision have an amplitude 
which is estimated to be of order Q.Gwh 2 ~ 10~ n at / = 3 x 10 _4 Hz, which is within reach of 
the planned sensitivity of LISA. For another range of parameters, we find that the amplitude is 
estimated to be of order Q,Gwh 2 ~ 10 -25 around / = 10 3 Hz, which is within reach of LIGO. Hence, 
it is possible to detect gravity waves from such a phase transition at two different detectors, with 
completely different amplitude and frequency ranges. 



I. INTRODUCTION 

Recently, the authors in pQ have analyzed the scalar 
field theory of a standard model Higgs field whose po- 
tential includes non-renormalizable operators up to mass 
dimension 8. In this model, a second true vacuum ap- 
pears which may be responsible for the observed late- 
time acceleration of the universe. While the electroweak 
phase transition remains second order, the subsequent 
phase transition from the second false vacuum to the 
true vacuum is first order, giving rise to the possibility 
of bubble nucleation. In this picture, we are sitting in 
the second false vacuum (located e.g. around <j> ~ 0.3 
TeV in Fig.[T]), awaiting the plunge into the true vacuum 
(around 4> ps 0.8 TeV). There are certain fore-bearers of 
this impending doom that are in principle observable. 

In particular, it is well-known that bubble collisions 
lead to the production of gravitational waves (GW), 
due to the breaking of symmetry. The strength of the 
phase transition, as well as the amplitude of the GW 
are strongly model and parameter dependent. Thus 
it seems worthwhile to investigate whether this non- 
renormalizable Higgs model can have a background of 
interest for current and proposed GW detectors 1 , specif- 
ically LIGO, Advanced LIGO, Geo 600, Virgo [2] and as 
well as the proposed LISA experiment [3]. 



The two main sources of GW production are the nu- 
cleation and collision of bubbles of true vacuum and the 
onset of turbulence. In this paper we shall restrict our- 
selves to the first case, that of collisions of bubbles of true 
vacuum. In order to tackle this problem, we rely on the 
pioneering work of [4, 5 who developed a semi-analytical 
model for the GW spectrum emitted from bubble colli- 
sions. 

This paper is organized as follows. Sectio n [TT] briefly 
reviews the model proposed in pQ. In Section |HI| we re- 
view the thin-wall approximation of 4. 5 . In Section IV 
we present the results of our calculation. Finally, we con- 
clude in Section fVl 



II. ELECTROWEAK MODEL 

In [1] the authors considered the standard Higgs po- 
tential with additional higher order non-renormalizable 
terms. These non-renormalizable terms do not pose a 
problem as long as they are suppressed by small enough 
coupling constants, i.e. large enough masses. The poten- 
tial considered is given by, 
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1 It should be pointed out that a possible detection of gravitational 
waves from bubble collisions of this phase transition would give 
us only a brief time before the joint bubble of true vacuum passes 
through us - the bubble wall moves at close to and gravitational 
waves move exactly at the speed of light. 



For later use, we rewrite the potential as 
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In the latter form, the finite temperature effective po- 
tential formulation of the theory is more easily written 
down. 

The parameter eg assumes the role of a controllable 
fine tuning of the potential. When £0 = 0, the potential 
in Eq. ^ has four degenerate minima at <f> = ±0i,±02- 
When £0 7^ 0, the degeneracy of the vacua between 4> = 
<Pi^<p2 {4> — ^0ii~</ , 2) is broken and there is an energy 
difference between the energy densities of the two vacua, 
as shown in Fig. [TJ (at T = 0). The difference in energy 
is given as 



^ = £0^(^-0?), 
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where <f>2 > 4>i- 

Finite-temperature effects are approximated by adding 
a thermal mass term to the potential in Eq. ([2]), hence 
the potential takes the form V(<j), T) = cT 2 cf> 2 /2+V(<f>, 0), 
where c is generated by the quadratic terms that acquire 
a 0-dependent mass in the high-temperature expansion of 
the one-loop thermal potential. In general there are also 
terms in the high-temperature thermal expansion that 
are proportional to T 4 c/> 4 . However, it was argued in [5] 
that these terms only lead to small corrections to the 
potential, therefore we shall ignore these contributions. 
In terms of the temperature, the effective potential then 
becomes 
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where 
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Using the procedure outlined in [71 H] , the constant c is 
found to be given by 
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where g and g' are the SU(2)l and U{1)y gauge cou- 
plings, and y t is the top Yukawa coupling. Using the 
excepted values of these couplings and that the Higgs 
mass is 0.08 TeV < m# < 0.15 TeV, one finds that from 
Eq. (10), c w 1.7. The critical temperature is defined 
when MT) = is zero or in other words 
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FIG. 1: The temperature dependence of the <j> Higgs field 
potential given in Eq. (JsJ) . The values used in this plot are 
4>i = 0.246 TeV, (j} 2 = 0.8 TeV, A 3 = 0.154 TeV -4 , Vq = 
and e = 0.01 TeV" 4 . 



In this case, the vacua are degenerate, as discussed above. 

In Fig. [T] we plot the temperature dependence of the 
potential in Eq. ([8]). Since the potential is symmetric 
for <j) — > —cj>, we will consider only the <j) > half-plane. 
For very high temperatures (T » T c ), the quadratic 
term of the potential dominates and the potential con- 
sists of a characteristic "U" shape with a single minimum 
at <j> = 0. Here, the Higgs field has zero expectation 
value and the electroweak symmetry is left unbroken. As 
T decreases, (j> = becomes a maximum (T > T c ) and 
the field begins to roll down the potential toward the 
first minimum, which is forming at <fi = cj)\. The second 
minimum is forming at <f> = <f>2 at higher energy, making 
it non-accessible since it is a less favorable state. The 
non-zero expectation value of the Higgs field at <p = <pi 
breaks the electroweak symmetry. This is the standard 
picture of electroweak symmetry breaking from a field 
going through a second order phase transition. However, 
the Higgs potential keeps evolving as the temperature 
decreases. The the second minimum at <j) — 2 drops 
to lower energies until it becomes the true global min- 
imum after the temperature falls below T c , causing the 
minimum at <j) = <fci to be the false minimum. In other 
words, the Higgs field sitting there in the false vacuum 
can undergo a first order phase transition by tunneling 
to the global minimum at 4> — <f>2 and forming bubbles 
of true vacuum. 

As is well known, both nucleation of bubbles and their 
collisions can source gravitational waves. In the next sec- 
tion we will review the production of gravitational waves 
via bubble collision. 



III. GRAVITATIONAL WAVE PRODUCTION 
IN PHASE TRANSITIONS 



We are interested in production of gravitational waves 
from a first-order phase transition of the Higgs field. The 
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Lagrangian for the Higgs field (p is given by 



bubbles is t = 0, [TU] showed that Eq. (15 1 leads to 
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where the potential possesses at least two non-degenerate 
local minima, such as that in Eq. ^ (as shown in Fig. [Th. 

In the following, we use a metric with signature (H — ). 

Classically, the false vacuum state is stable: a field sit- 
ting in this vacuum will remain there forever. However, 
quantum effects can cause the false vacuum to decay to 
the true vacuum. This decay proceeds via the nucle- 
ation and expansion of bubbles of true vacuum, which 
spontaneously appear in the false vacuum. The bubbles 
then expand due to the energy difference between the two 
vacua, which induces an effective pressure on the bubble 
wall. This causes a bubble to expand with an initial ac- 
celeration, with its expansion speed quickly approaching 
the speed of light. 

In [S] , Coleman showed that the bubble evolves accord- 
ing to a Klein-Gordon equation with 0(3, 1) symmetry 
which implies that the position of the bubble wall is given 
by 
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where i?o is the initial radius of the bubble and x wa u 
denotes a fiducial point with in the bubble wall. When t 
is sufficiently large, we see that the position of the bubble 
wall is given by R(t) w t, hence the radius of the bubble 



wall is proportional to the elapsed time. From Eq. (12 1, 
the stress-energy tensor associated with the bubble is 



T M „(a?, t) = d^dyip - g^L 
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where the energy density of the bubble is the Too com- 
ponent. 

To determine the gravity-wave production we shall em- 
ploy the envelope approximation, originally developed by 
[lOj in the context of the linearized gravity approxima- 
tion. This approximation is valid for bubble sizes less 
than H . The energy radiated in gravitational waves is 
given by the Fourier transform of the spatial components 
of the stress-energy tensor in Eq. (14). Here, we will 



neglect the g^L term, since it is a pure trace term and 
does not contribute to the production of the gravitational 
waves. Following Weinberg (1 1 j . the Fourier transform 
can be shown to be 

T l3 {k,u) = ^- j dt e iut j datditpdjtp e-^H (15) 



where k is a unit wave vector. In the envelope approxima- 
tion, one assumes that the regions that have overlapped 
do not contribute substantially to the gravitational radi- 
ation. Thus, these regions are excluded from the spatial 
integrals. Assuming that the time of nucleation for the 
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where the sum is over the number of colliding bubbles. 
From Weinberg, the total energy radiated in gravity 
waves is given as 
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where Aijj m is the projection tensor for gravity waves: 

A-ij.im{k) = SnSj m 2kjk m Sn + —kikjkik m 

- 2%^m + ^Sijkik m + -Simkikj. (18) 

To compare our results with those found in [10] . we 
compute the case of two-bubbles colliding in the lin- 
earized gravity approximation. In Figs. [2] and [3] we plot 
some of our numerical results. Here d is the initial dis- 
tance between the two-bubbles. To simplify notation, we 
are plotting in units of Gp vac /2> which is just a constant, 
thus the relative magnitude of the curves in Figs. [2] and 
[3] will be reduced. In Fig. [2] we plot the scaled energy 
spectrum for initial separation of d = 60, while in Fig. [3] 
we plot the scaled energy spectrum per octave frequency 
interval. Comparing our results, we see that the results 
are qualitatively the same, hence consistent with those 
in [TO]. This is expected since we are using the envelope 
method, where the only dependence on the specific mod- 
els potential comes in the form of p vac , which is just a 
constant. 
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FIG. 2: Scaled energy spectrum for initial bubble separation 
of d = 60. 
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FIG. 3: Scaled energy spectrum per octave frequency interval. 



IV. GRAVITATIONAL WAVE PARAMETERS 

To estimate the spectrum of gravitational waves that 
are sourced by this phase transition, we need to estimate 
the tunneling rate for the decay of the false vacuum state 
to the true vacuum state as well as the temperature scale. 
We also need to go into some details of the subsequent 
evolution of the growth of the bubbles. In particular, we 
need an estimate of the velocity of the detonation front 

c 



A. Tunneling Rate and Temperature 

The decay rate is suppressed by the exponential of the 
effective action, T = TQe~ SB ^\ The time scale for the 
decay /3 given by 



dS E 



dt 



t=t, 



where i* is the time when the phase transition occurs. 
This gives 
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where t = it is the Euclidean (or Wick rotated) time 
coordinate. The full, 0(4) symmetric solution of Eq. ( plj ) 
and the nucleation of bubbles was analyzed first by [9"Ill2] 
for a universe at zero temperature. At finite temperature 
T, [E] pointed out that the field theory should be taken 
periodic in r with period T -1 . Hence Eq. (21) must be 
written as 
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where V(ip, T) is the temperature dependent effective po- 
tential. At sufficiently large temperatures, i.e. when the 
integrand is approximately time independent over time 
scales T _1 , the integration over r is reduced simply to 
the multiplication of T _1 , or S E = 5 3 /T 14J. Here 5 3 
has 0(3) symmetry and is given by 
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(jjj where the last equality assumes the thin-wall approxima- 



tion. Si (T) is the solution of Eq. ( 22 ) to zeroth order in 
«^(T) = E (T>?^I-^) 
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which is independent of the temperature as all T depen- 
dence is in SV. 

The temperature dependent bubble radius 
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is obtained by minimizing S3 . Plugging this back in ( 24 ) 
we find 
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Note that for this thin wall approximation to hold, the 
bubble radius must be much larger than the thickness of 



(19) the wall, R > (8$y(ip)) 



We can now estimate 
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16^ 
3 

16ir 2 



Sfd T 
Sf 



1 



§V(T) 2 T 
AcT 2 ip 2 - 8V{T) 
T 2 ip 4 5V(T) 



(29) 



The temperature at which the transition takes place 
is computed in a two step process. 

First, we fix eo by requiring that the probability is 
much smaller than unity for our Hubble patch to already 
have performed the second phase transition during the 
lifetime of the universe or in other words, by demanding 
that our Hubble patch is still in the false vacuum at ip\. 
Therefore we have 
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where tp is the Planck time. The smallest value of the 
action, i.e. the highest probability for the transition from 
4>i — > 4>2 to occur, is realized at zero temperature. Thus 
we can substitute the action at temperature T by the 
zero temperature action to obtain a lower limit on the 
action. 
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which should be larger than the critical action 
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Then, we estimate the transition temperature by 
equating 
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and solving for T, which should at least give us a rough 
estimate. 

Taking tp 1 = 0.246TeV, <p 2 = 0.8TeV, A 3 = 
0.154TeV~ 4 and S ait = 425, we find e = 0.15TeV" 4 
and a transition temperature of order = 0.0128TeV. 

From the transition temperature we can now determine 
/3/H*, the quantity of interest (as we shall see in the next 
section). We first need to make use of the fact that during 
radiation domination 
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where N(T) is the effective number of particles, see e.g. 
15J for more details. Thus we have 
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B. Growth of Bubbles 

In this subsection, we estimate several quantities per- 
taining to the growth of the bubbles, notably the rela- 
tivistic 7 of the bubble wall. 

An expanding bubble of true vacuum behaves very 
much like a detonation. This was first examined in detail 
by 02] and the to first order phase transitions by [S]. 



In a detonation, high temperature gas enters the nu- 
cleated bubble with supersonic velocity. This prevents a 
shock from preceding the bubble wall. The velocity of 
the detonation front is given by 
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where the relativistic 7 is approximately the ratio of vac- 
uum energy to thermal energy. There are several ways 
to estimate its value, but they all amount to numbers of 
the same magnitude. 

The efficiency factor, or fraction of latent heat that 
goes into the kinetic energy rather than thermal energy, 
is defined as 



K{i) = 
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(41) 



To determine the latent heat and vacuum energy asso- 
ciated with the transition, we begin with the value of the 
potential at the broken phase minimum which is also the 
difference in free energy density between the two states 
of the system. We denote this by 



B{T) 
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where v(T) is the value of the broken phase minimum. 
The vacuum energy associated with the transition is then 



£* = BIT)- TB'(T) 

T=T, 

This results in a relativistic factor 7 of 
30£* 
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be 



Alternatively, we could estimate 7 just like in |10j to 
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(45) 



Both estimates for 7 resultin a very large number. The 
wall velocity then approaches speed of light, and the ef- 
ficiency approaches unity 



(46) 



Hence we are justified in using the approach in Sec. |III| 
One final remark is in order. In Sec. [TT] we have ar- 
bitrarily set Vq = 0. However, the chosen value of Vq 
effects the value of the vacuum energy associated with 
the transition, see Eq. (43). In particular, the presence 



of a non-zero Vq acts as a overall shift in the vacuum en- 
ergy, since the derivative of the potential at the broken 
minimum, B'(T), is independent of Vq. Therefore one 
would write the vacuum energy as £ * = + Vq 1 . How- 
ever, since 7 is so large, the constant shift, which is on 
the order of a few TeV, would increase the above num- 
bers. But as 7, £d,i and K are already indistinguishable 
from their values in the limit 7 — > 00 , this shift is not im- 
portant for the determination of the gravitational wave 
spectrum. 
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C. Gravitational Wave Spectrum 

The spectrum of gravitational waves emitted from bub- 
ble collisions was first found by [5] . The fraction of the 
energy density compared to the critical energy density is 
given by 
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and the characteristic amplitude of a gravity wave is ap- 
proximately 
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In the previous section we found that 7 — > 00 so that 
£d — » 1 and k — > 1. Plugging in these values, we find 
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(52) 



Therefore for bubble collisions we get O/i 2 ps 1.7 x 10 11 
and /i c ps 3.1 x 10~ 19 , peaking at a frequency around 

fmax « 3.5 X lO" 4 HZ. 

We briefly note that for smaller values of £0 ~ 
0.005TeV~ 4 , the transition occurs at temperature ps 
3GeV. This gives /3/-H* ps 3.7 x 10 9 , leading to an am- 
plitude r^cw/i 2 ~ 1.1 x 10 -25 at peak frequency f max ~ 
10 3 Hz. 



V. CONCLUSION 

We have considered the standard Higgs field, 
albeit with a potential containing higher order 
non-renormalizable terms. Introducing these non- 
renormalizable terms induced a second minimum in the 
potential, see Fig. [l] The presence of the first quadratic 



term in Eq. (JS|) ensures that the two minima are non- 
degenerate at sufficiently low temperatures. There is a 
false-vacuum (a local minimum) and a true-vacuum (a 
global minimum) . In the context of this model we assume 
that our universe is currently in the false vacuum. The 
presence of the second lower minimum means that the 
false- vacuum can decay via the nucleation and expansion 
of bubbles of true- vacuum from a first-order phase tran- 
sition. As the bubbles are growing, they will eventually 
collide with other bubbles, resulting in the production 
of gravity waves. In this paper we have investigated the 
gravity wave spectrum from bubble collisions using the 
envelope approximation developed in [5]- We have also 
considered the resulting gravity waves from detonating 
bubbles. 

For bubble collision from a given first-order phase tran- 
sition, knowledge of the model-specific parameters (3, £, 
k and 7 suffices to determine the resulting gravity wave 
spectrum, see Eqs. (47) - (49). The effective potential 



for bubble nucleation suffices to determine the spectrum 
of gravitational radiation. Here we find that for the 
Higgs potential containing the non-renormalizable terms 
Eq. ([2]), the transition occurs at a temperature ps 
12.8GeV (see Section |IVA |), » 300, and 7 - 10 4 , 

so that £d ~ 1 fa k (see Section IV B ) . Using these results 
and Eqs. p7b - |49) from Section |IVC] we find that the 
energy density of gravity waves Slh z ps 1.7 x 10~ n , their 
characteristic amplitude h c ss 3.1 x 10~ 19 , and a peak 
frequency around f m ax ~ 3.5 x 10~ 4 Hz. 

The question then arises whether or not the gravity 
wave signal is potentially detectable by current or fu- 
ture gravity wave experiments. Here we consider some 
of the possible detectors and their projected sensitivities 
and compare with our results. First we shall consider a 
current land-based detector. From [T51 Q15] the ultimate 
sensitivity to a background for the Laser Interferomctric 
Gravitational Wave Observatory (LIGO) is an amplitude 
around 2 x 10~ 25 at a frequency of 100Hz. However, as we 
saw in Section IV C we found that f max ps 8.2 x 10 _4 Hz 
with an amplitude of 1.7xl0 -11 , which are well below 
this requirement. Therefore a detection by LIGO is very 
unlikely. Next we consider a future space-based exper- 
iment. The peak sensitivity for the Laser Interferome- 
ter Space Antenna (LISA) [16l [20l [21] is reached in the 
frequency range between 3 x 10" 2 Hz and 10 _1 Hz with 

10- 10 
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From Section IV C we see that 



our results fall within this range of detection from LISA, 
hence a detection by LISA could be possible given our 
results. 

One thing to note, however, is that the value of £0 = 
0.15 TeV -4 is an upper limit since smaller values of Eq 
further decrease the bubble nucleation rate. Repeat- 
ing the calculations above, we see from Eqs. ( |47| - ( [49] ) 
that as £0 decreases, Vlcwh 2 and h c (f max ) both decrease 
while f m ax increases. For the value of £0 given in Fig. [T] 
&Gwh 2 ~ 3 x 10~ 18 which is well outside the range of 
LISA. However, as £0 falls below 0.005TeV^ 4 the peak 
frequency and amplitude shift into a window that may 
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be detectable by LIGO. For this value of Eq the transi- 
tion occurs at temperature T* 3GeV, /3/H* w 3.7x 10 9 . 
This lead to an amplitude of G( y/! 2 ~ 1.1 x 1CU 25 at 
peak frequency f max ~ 10 3 Hz. 

Even though we don't consider these here, additional 
contributions to the gravity waves spectrum (which could 
be comparable to, or maybe even stronger than bubble 
collisions) could come from the period of turbulence after 
the phase transition, see [52]. First of all, turbulence it- 
self sources gravitational waves. In addition to that, tur- 
bulence can also source magnetic fields during the time 
of the phase transition, see [551 H3] • These induced mag- 
netic field will then act as a source of gravitational ra- 
diation until they are damped out. The gravitational 
radiation from the latter process is subdominant to that 
of the turbulence, however, its peak frequency is much 



higher and could potentially be detected at a land-based 
detector such as LIGO. It would be interesting to inves- 
tigate the spectrum of gravity waves sourced by these 
effects as well as the induced magnetic fields. 
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